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Abstract 

We find that single spin asymmetries in NRQCD are non-vanishing 
in general. They are proportional to the imaginary parts of some non- 
perturbative matrix elements. With statistics of about 10 6 identified 
J/ijj's, or 10 5 identified X2, it is possible to measure these imaginary 
parts even if they are an order of magnitude smaller than the real 
parts. Such statistics are quite reasonable at polarised HERA N and 
other future experiments. 
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1 Introduction 



Single polarised initial states can give rise to asymmetries in two ways. If 
a final state polarisation is observed, then the difference between the cross 
sections when the initial and the final state spins are parallel and anti-parallel 
gives an asymmetry. These do not yield information beyond what comes from 
double polarised initial states. Recent interest has focussed on single spin 
asymmetries in which final state spins are not observed — 

s = w = ^Mfc) 

where Vex is the difference of cross sections with the initial particle polarised 
in opposite directions. Such asymmetries do not exist in leading twist compu- 
tations in perturbative QCD. They seem to arise either from next-to-leading 
twist-3 effects, or through non-perturbative effects like intrinsic transverse 
momenta of partons. 

A formalism for computing the production cross-sections of heavy quarko- 
nia, called non-relativistic QCD (NRQCD) |]], has been developed. In 
NRQCD the non-perturbative long-distance physics is factored from the short 
distance perturbative physics. This seems like an ideal place to search for 



single polarised asymmetries of the kind in eq. (ET]). In this paper we re- 
port a set of computations which shows that such asymmetries indeed exist 
in NRQCD, yield new and interesting non-perturbative information about 
quarkonia, and are measurable. 

NRQCD is a low-energy effective theory for quarkonia. The action is 
written in terms of all possible operators consistent with the symmetries of 
QCD. All momenta in NRQCD are cut off by some scale A. The coupling 
associated with each term is identified through a perturbative matching pro- 
cedure. The applicability of NRQCD to quarkonium production depends on 
the proof that final state effects factorise in the hadronisation of a heavy 
quark pair to a quarkonium. Such a proof has been given for production at 
large transverse momenta jl] and successful phenomenology has been done 
[0]. For low-energy production of quarkonia a proof is lacking, but a rea- 
sonable phenomenology arises when cross sections are computed assuming 
factorisation 0, |], Double polarisation asymmetries have been computed 
for quarkonia in this framework || [7|, |J. Numerical results for single polari- 
sation asymmetries computed in the colour singlet model have been quoted 
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in the literature ||, and preliminary results in NRQCD obtained in ||. 

The NRQCD factorisation formula for the inclusive production of heavy 
quarkonium resonances H with 4-momentum P can be written as 

d -= ^(^^^^(W^J)' (L2) 



rfV(T= - 7 —^—^ ij C ij (jC i n(H)K]). (1.3) 



1 d 3 P 
$ (2n) 3 2E r 

where $ is a flux factor. The coefficient functions and Cy are computable 
in perturbative QCD and hence have an expansion in the strong coupling a s 
(evaluated at the NRQCD cutoff). Although each matrix element in the 
sum above is non-perturbative, it has a fixed scaling dimension in the quark 
velocity v. Then the NRQCD cross section is a double series in a s and v 2 . 
For charmonium a numerical coincidence, a s ~ v 2 , makes this expansion a 
delicate one to handle. For bottomonium, the series is easier to analyse since 
v 2 < a s . 

The fermion bilinear operators /Q are built out of heavy quark fields 
sandwiching colour and spin matrices and the covariant derivative D. The 
composite labels i and j include the colour index a, the spin quantum number 
S, the number of derivative operators N, the orbital angular momentum L, 
the total angular momentum J and the helicity J z . The hadron projection 
operator 

11(H) = Y,\H,s)(H,s\, (1.4) 

s 

(where s denotes hadron states with energy less than the NRQCD cutoff), 
is diagonal in L, S, J and, when the final state spin is unobserved, connects 
states of opposite J z . As a result, the operators K{ and K.j in eq. ( |1.3j ) are 
restricted to have equal L, S, J and opposite J z , for both da and dVa. 

The ^-dependence of these matrix elements can be factored out using 
the Wigner-Eckart theorem — 

(zmmh = ^jO»e s+i L»), (1.5) 
(wmict) + h .c. = ^r^^n (i-6) 

The factors of 1/(2 J + 1) come from a Clebsch-Gordan coefficient and are 
conventionally included in the coefficient function. Each diagonal matrix 
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element, O, is the product of a reduced matrix element (if||/Cj||0) with its 
complex conjugate (0| \H), and is real. However, the off-diagonal matrix 
elements, V, may have non-vanishing phase. We introduce the notation 

(m(H)JC]) - h.c. = ^^C S+l L N j, 2S+1 Lf), (f .7) 

where the Wigner-Eckart theorem has been used as before to factor out the 
J z dependence of these matrix elements. The NRQCD power counting rule 
for all three types of matrix elements is — 

d = 3 + N + N' + 2(E d + 2M d ), (1.8) 

where Ed and Md are the number of colour electric and magnetic transitions 
required to connect the hadronic state to the state /Q|0). 

In section ^| of this paper we present the main results for Vex after a 
telegraphic review of the threshold expansion technique ||10|| , the kinematics 
and the appropriate Taylor series expansion of the perturbative matrix el- 
ement ||11|| . Although the results presented in this paper are for quarkonia 
produced at small transverse momenta, it should be clear that similar effects 
should arise also at large transverse momenta. Observable consequences are 
discussed in the final section ||. 



2 Computing the Asymmetry 

We choose to construct the coefficient functions using the "threshold expan- 



sion" technique of [10|. This consists of calculating, in perturbative QCD, 
the matrix element M. connecting the initial states to final states with a 
heavy quark- ant iquark pair (QQ), and Taylor expanding the result in the 
relative momentum of the pair, q, after performing a non-relativistic reduc- 
tion of the Dirac spinors. The resulting expression is squared and matched to 
the NRQCD formula of eq. (|1.3|) by inserting a perturbative projector onto 
a non-relativistic QQ state between the two spinor bilinears. The coefficient 
of this matrix element is the required coefficient function. 

In this paper we evaluate the polarised cross sections to order a 2 s v 9 . This 
requires a Taylor expansion to order, N + N' < 6, as can be seen by setting 
d = 9 and Ed = Md = in eq. ( |1.8|) . A simplification occurs because the 
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perturbative projector has only one term — 

U(QQ) = \QQ)(QQ\. (2.1) 

In agreement with || [10], [11] we use the relativistic normalisation of states 

(Q(p,Z)Q(q,v)\Q(p , ,e)Q(q , ,v')) = ±E p E q (2Ttf5\ V - V ')5\q-q'), (2.2) 

with the spinor normalisations = rfrj = 1. Expanding E p = E q = 
y/m 2 + q 2 in q 2 allows us to write the spinor bilinears in terms of transition 
operators built out of the heavy quark field. 

The kinematics is very simple to leading order in a s . The momenta of 
the initial particles are pi and p 2 . We take pi to correspond to the polarised 
particle and assume that it lies in the positive z-direction and that p 2 is 
oppositely directed. The net momentum P = pi + p 2 . The 4-momenta of Q 
and Q (p and p respectively) are written as 

p = ^P + Ljq 1 and p = l -P - Lrf '. (2.3) 

Note that p 2 = p 2 = m 2 , where m is the mass of the heavy quark. The 
space-like vector q is defined in the rest frame of the pair, and Lj boosts it to 
any frame. We shall use Greek indices for Lorentz tensors and Latin indices 
for Euclidean 3-tensors. A property of the boost matrix that we shall use 
many times is 

M 2 

^xpPiP^Li = —eijkZk, (2.4) 

where z is the unit vector in the ^-direction. This can be derived from some 
of the identities listed in |H| and the kinematics given here. 



a ±1 = T— (a x ± iay) , = a z , (2.5) 



The power counting rule in eq. (|1.8|) requires that we express the Fermion 
bilinear operators as spherical tensors. Recall that any vector can be 
written as a spherical tensor of rank 1, with the components 

1 

V2 

where the subscripts ±1, are helicity indices. In NRQCD higher rank 
spherical tensors are constructed by coupling such rank 1 tensors successively 
[PH| . Some useful identities are 

djbj = a b - (a+i&_i + a_i& + i) = — \/3[a, (2.6) 
iej k iajb k zi = (a+i6_i - a_i6 +1 ). (2.7) 
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We have introduced the notation [a, h] J M to denote two spherical tensors a and 
b coupled to total rank J and helicity M. The coefficient of the term [a, 6]§ in 
eq. ( |2.(j| ) can be obtained from the appropriate Clebsch-Gordan coefficients. 



2.1 QQ 

As a simple example of the techniques used, we write down the matrix ele- 
ment for the subprocess qq — > QQ — 



• 2 

M = -j^iv(p 2 ,) ltl T a u( Pl ,h)}L» 

x 



4 

M^a 3 T a 7i q^Hq • a)T a r/ 



(2.8) 



where T a is a colour generator, u and u are the light quark spinors and £ 
and 77 are the heavy quark Pauli spinors. The equations of motion for the 
initial state quarks have been used to obtain the explicitly gauge invariant 
matrix element above. The desired Taylor series expansion is written down 
using the relation M 2 = 4(m 2 + q 2 ) to expand all factors with M. 

The squared matrix element for this process is trivial. The light quark 
spinor factors require the projection operator + /ry 5 )/2. The polarised 
cross section is obtained from the part proportional to h. Doing the Dirac 
algebra and using the identity in eq. fl2.4p to simplify the result, we find— 



• 2 

V|-M| a = -^I^kikMjMl (2.9) 

where Mj is the heavy-quark spinor bilinear in eq. ( |2.8|) . Using eq. (|2.7Q 
and the Wigner-Eckart theorem, it is easy to see that the contribution of the 
diagonal operators vanishes. 

The final results for the cross section differences, to order v 9 , are — 



Jli> 



v u qq 



V(T- C 
<?<2 

54m 6 
54m 6 



5(s — Am 2 
5(s — 4m 2 



>7f l 

2 

— -T XJ [ 
V3 8 



CO 3q2' 



3c0 3c2- 



7V5_ 
12m 2 



(2.10) 
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where s is the parton CM energy. The expressions are exactly the same 
whether the initial quark or anti-quark is polarised. The matrix elements 
are all of order v 9 . The corresponding unpolarised cross sections have been 



written down in |TT|. They are of lower order in v. 

The fact that the asymmetry is non-zero for the 3 Po state, Xo, might come 
as a surprise. However, one should note that the matrix element involved 
in its production entails the radiation of at least one soft gluon. Thus, in 
qq annihilation a xo can only be produced along with some light hadrons. 
The total angular momentum of the final state is 1, making it possible to 
generate this single spin asymmetry. 

2.2 gg QQ 

The squared matrix element for the gg process is technically a little more 
complicated. We work in a class of ghost-free gauges called planar gauges 
|H| . The density matrix for an initial state gluon of momentum p and helicity 



h in these gauges is given by 



e ^(p, h)et(p, h) 



p ■ V p ■ V 



(2.11) 

The vector V defines the gauge choice. We write V = C\P\ + c 2 p 2 , with 
ci/ca ~ 0(1). 

The sum of the matrix elements arising from the three Feynman diagrams 
(s-channel gluon exchange, A4 S , and t- and ^-channel quark exchanges, Ait 
and M. u ) can be decomposed into three colour amplitudes— 

M = \g 2 5 ab S + l -g 2 d abc D c + l -g 2 f abc F c . (2.12) 

The colour amplitudes S and D involve only Ait + -M-u, whereas F involves 
Ai s as well as Ai t — Ai u . 

In order to write down our results, we find it convenient to introduce the 
notation 

A = j^£\ai,uPiP2 e i e 2 and s ij = A i z j +A j Zi-B i j + e 1 -e 2 z i z j , 

(2.13) 
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where 

A = h (ei • Li€ 2 ■ pi - e 2 • Uei ■ p 2 ) , 

(2.14) 

Bij = €\ ■ Li€2 ■ Lj + 62 ■ Li€\ ■ Lj. 

Here ej is the polarisation vector for the initial gluon of momentum p^. In 
order to identify all terms to order v 9 we need the colour amplitude S to 
order q 5 — 

S = - (8jm) A ^V) + ifS rm (q^v) ~ (^) Az m z n (q m q^v) 

(q m q n q p ^(jiri) 



, 16 



64 



^jm^n^y -j- 2/77. 3 m n V 

Az m z n z p z r (q m q n qPq r ^r ] ) 
M 



S jm z n z p - ^ + 2m ^- m 5 np j z r z s (q m q n q p q r q s ^a J r]). 

(2.15) 

The amplitude -D differs only through having colour octet matrix elements 
in place of the colour singlet ones shown above. For the colour amplitude F 
we need the expansion 

F c = ~ (^?) Az m (q" 1 ^) + ^mZn (q m q n e^T c V ). (2.16) 

In all three colour amplitudes, the terms in A are spin singlet and those in 
S are spin triplet. 

The density matrix in eq. ( |2.11| ) yields 

VS-S* = ^-^ ha i h rnS jm S], m ,cyd ml = { [a, b] 2 _ 2 [c, d)j - [a, b)\ [c, df_ 2 } , 
4 h 

(2.17) 

where a, b, c and d are Euclidean 3- vectors, written in terms of the Euclidean 
components on the left and as spherical tensors on the right. Gauge invari- 
ance is obvious from the fact that the final result does not depend on the 
gauge choice V. In addition, we find the other gauge invariant contractions, 
VA ■ A* = and V.4. ■ S* ^ 0. Since A always comes with a spin-singlet 
operator and S with a spin-triplet, this last combination always gives rise to 
a product of two fermion bilinears of opposite parity. Such operators vanish 
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in NRQCD. Finally, the single spin asymmetries can be found by retaining 
only the terms in S in the three colour amplitudes S, D and F. 

After retaining only those spherical tensors which contribute to order v 9 
(they are tabulated in |fLl|| ) and dropping all the diagonal tensors, very few- 
terms remain. It is clear, for example, that the colour amplitude F could 
contribute through only one term. However, its contribution must then be 
diagonal, and as a result this colour amplitude cannot give a single spin 
asymmetry at this order. Similar arguments tell us that only two terms 
contribute from each of the S and D colour amplitudes. The coefficient 
functions can be read off from earlier computations of the unpolarised cross 
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The direct J ftp single polarised cross section difference is 
Va J g ^(s) = ^0^(9), where <p = ^6(s - 4m 2 ), (2.18) 



and & J J^ (d) denotes combinations of non-perturbative matrix elements from 
the colour amplitude a (= S, D or F) at order v d for the cross section 
difference. The matrix element required here is 

9^(9) = - *Tt/*('Pl*i#). (2.19) 



The unpolarised cross section [11| starts at order v 7 . The cross section dif- 
ferences for ip', T and all other 3 Si states are given by the same formulae; 
only the appropriate matrix elements have to be plugged into eq. ( |2.19|) . 



Upto order v 9 , the single polarised subprocess cross section differences for 
1 5'o, l P\ and 3 Po,i quarkonia vanish — 

Va*°(s) = Vaf g {s) = V^(S) = V6* (*) = 0. (2.20) 

The 3 P2 cross section difference is 

V^(5) = <p (1) [eT(7) + Qf(9)} (2.21) 

where the combinations of non-perturbative matrix elements are 



(2.22) 

The unpolarised cross section JTTj] starts at order v 5 . 
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2.3 ig -> QQ and 77 -»■ QQ 

The matrix elements for the two processes jg — > QQ and 77 — > are 
closely related to the gg amplitudes. It is obvious that 

M ig = geD, and = e 2 S, (2.23) 



where D and S are the colour amplitudes given in eq. (|2.12|) , and e is the 
charge of the heavy quark. 

The jg cross sections for the production of any quarkonium state can be 
obtained from those for the gg process, by the prescription — replace a 2 s by 
aa s , delete the 65- and 9^ terms, and replace the colour factor 5/48 for the 
terms in by 2. The 77 cross sections are obtained with the prescription — 
replace a 2 s in by a 2 , delete the 0£> and 9^ terms, and replace the colour 
factor 1/18 for the terms in G5 by 16. Only 3 Si quarkonia have non- vanishing 
single spin asymmetries in almost-elastic photo-production, and only \2 can 
have non-zero asymmetry in photon-photon fusion. 



3 Experimental Outlook 

The observable asymmetries in hadro-production are obtained by polarising 
one of the initial hadrons. These quantities are obtained by convoluting the 
cross sections found in Section |] with appropriate parton density functions. 
It is easy to see that 

s h {s s = Vfrg(a)V£ft(s, y) + Vag(g)V£ flfl (s, y) 

where y/s is the centre of mass (CM) energy of the colliding protons and y is 
the rapidity at which the final state quarkonium, H, is observed. With the 
convention that the positive z-direction is defined by the direction of motion 
of the polarised initial state hadron, the parton luminosities can be written 

as 

V A<? = J2 A( lf {xi)q f {x 2 ) + Aq f (x 1 )q f (x 2 ), VC gg = Ag(x 1 )g(x 2 ), 

(3.2) 
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Figure 1: The ratios C qq /C gg (full line), VC qq /C gg (dashed line) and 
VC gg /Cg g (dotted line) computed using the GRSV LO parton density set 
function of the CM energy a/s at zero rapidity. 



where x\ = {2m/ ^/s) exp(y) and £2 = (2m/-y/i) exp(— y) are the fractional 
momenta of the polarised and unpolarised initial hadrons (respectively) car- 
ried by the partons. The quantities Agy, Aqj and Ag are the usual polarised 
quark, anti-quark (of flavour /) and gluon densities, and the corresponding 
symbols without the A stand for the unpolarised densities. The unpolarised 
parton luminosities C qq and C gg are given by similar expressions where the 
polarised parton densities are replaced by the unpolarised densities. 

Numerical values of the ratios of several of these luminosities are plotted 
in Figure [l] as a function of \fs. It is clear from the figure that C gg is the 
largest and that both the polarised and unpolarised qq luminosities are small. 

The asymmetries Spp b , and Spp b are likely to be too small to measure. 
The ratio of the non-perturbative matrix elements which determines the 
analysing power of these processes are of order v 6 . Since v 2 ~ 0.3 for char- 
monia and about 0.1 for bottomonia, these asymmetries are unlikely to be 
observable. For both and non-vanishing contributions at order v 9 
come only from the qq process whereas the cross sections are dominated by 
the gg process. Thus, these asymmetries are small because of the smallness 
of the ratio V C gq / C gg . is further suppressed by an analysing power which 
is of order v 4 . 
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The most promising observable is the single spin asymmetry for X2~ 



Sp p ( s. 



y) 



5/3J? 



(3 j> \ 3 



Pi) 



Ag(zi) 

gfa) 



(3.3) 



The analysing power is of order v 2 ; the order v 4 correction can be easily 
written down using eq. (|2.22|) , if required. Every matrix element can be 
written as a dimensionless number multiplying appropriate powers of v and 
the NRQCD cutoff A ~ m. For the real parts of the matrix elements, there 
is some evidence || that this dimensionless number is independent of which 
matrix element it comes from, and depends only on the quarkonium state. 
We use this assumption along with the definition, 



ire pi, 3 pi 



tan$ X2 Pf ( 3 P 2 \ 3 P 2 3 



(3.4) 



of the phase of this operator^ to ask what is the expected sensitivity of 
experimental measurements of $ X2 . 

First, we rewrite the asymmetry in the form 



Next we note that the measured asymmetry, 



5/3tan$ X2 Agfa) 



l + v\/5/3 9 fa) 



N+ + N. 



(3.5) 



(3.6) 



(iV + is the number of events observed when the initial hadron is polarised 
in the positive direction, and iV_ when it is negative) is smallest and most 
susceptible to error when N + ~ iV_ = N/2; and the expected error in S is 
less than or equal to 1/y/N (the equality is reached when the errors in iV + 



and iV_ are perfectly correlated). By inverting eq. (3J3) we can readily see 
that, near S = 0, the sensitivity of the experiment is 



|tan$ X2 | > 



1 



« 2 V^/3 9fa) 



N5/3 A 9( 



(3.7) 



"If the scaling rule in eq. (l.E) is to determine the importance of various operators, 



then we should expect that the phase is independent of the operator and depends only on 
the hadron. 
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Recent experiments at yfs 



38.8 GeV have observed a J/ip cross section 
About 15% of these come 



of about 400 nb in proton nucleon collisions JT4 
from decays of \2- At the future polarised-HERA N the CM energy is likely 
to be about 39.3 GeV with an integrated luminosity of 80 pb~ x per year [|]. 
Even if only the forward moving J/^'s are recorded, this should yield about 
a million observed J/ip in the dimuon channel per year of run, and hence 
about 10 5 identified X2, assuming a reconstruction efficiency of about 67% 
for the x's. The appropriate ratio of parton densities^ is around 0.4. If the 
degree of polarisation of the target is about 0.8, polarised-HERA N should 
be able to measure 



tan$ X2 | > 0.03 



(single polarised - HERA N). (3.8) 



Thus, single spin asymmetries in \2 should be measurable at polarised-HERA 
N even if the imaginary part of the matrix element is more than an order of 
magnitude smaller than the real part. Current Fermilab experiments have 
statistics of about 10 5 observed J/ip, and hence about 10000 X2 (assuming 
a reconstruction efficiency of 67%). With similar target polarisation their 
sensitivity is only a factor 3 less. 

For direct production of 3 Si quarkonia the single spin asymmetry is given 

by 

' Agfa) 



qlp'. 



y) 



s£(9) 



e£(7) 



eg(9) + (9/5)e£(9) 



9 fa) 



(3.9) 



where the unpolarised matrix elements are given in 
assumptions as before, we can rewrite this as 



*S , pp( s ) y) 



0.65^^11$,/,/ 
1 + l.lu 2 



Agfa 
gfa) 



. With the same 



(3.10) 



The analysis is very similar to that for \2- However, the results are much 
weaker, since the ip' cross sections are two orders of magnitude smaller. In 
this case the polarised-HERA N and Fermilab experiments can only measure 
an imaginary part which is no smaller than a third of the real part. 

The asymmetry for J/ip is the hardest to predict, due to feed down from 
radiative decays of x an d ip'- The problem can be simplified by working with 



3 By our convention, x\ is the momentum fraction of the parton from the polarised 
target, and xi from the unpolarised beam. 
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a data set consisting of directly produced J / ifj. Then the analysis follows from 
the analysis of ip' above. A million J/^'s can easily be obtained at polarised- 
HERA N. The sensitivity is similar to that for \2- At Fermilab fixed target 
experiments statistics are poorer by a factor of 10, and the sensitivity is 
smaller by a factor of 3. 

We conclude with a summary of our results. NRQCD predicts non- 
vanishing single polarised asymmetries for 3 Si and 3 P 2 quarkonia. The 
asymmetries are proportional to imaginary parts of some off-diagonal non- 
perturbative matrix elements. Since these are unknown, it is not possible to 
predict the values of these asymmetries. However, at polarised-HERA N and 
with polarised targets at Fermilab it is possible to measure the values of these 
imaginary parts even if they are an order of magnitude smaller than the real 
parts. Such experiments are extremely exciting, since they will probe an as 
yet unknown sector of the theory of quarkonia. 
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